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We present next-to-leading order QCD predictions for single-inclusive and di-jet observables, relevant to the 
collisions of two on-shell photons, obtained with a recently completed computation based upon the subtraction 
method. We compare these predictions with available theoretical results, with LEP data, and we discuss the 
various uncertainties which affect them. 



1. Introduction 

During the past two years different analyses 
|1I2I3I4| of jet and di-jet production from real 
photons scattering at LEP have been presented. 
A unique computation it is available in the liter- 
ature |5I6| at NLO for this process. This com- 
putation is based on the slicing method to can- 
cel analytically virtual and infrared 1 divergences. 
This technique relies on an approximate treat- 
ment of the phase space and of the matrix ele- 
ments, which induces cancellations between large 
numbers. Hence in particular regions of the phase 
space it is possible that these cancellations are 
difficult to achieve numerically. To avoid this 
problem (it could be substantial for specific exclu- 
sive observables which have been measured pQ), 
it seems appropriate to use an exact technique to 
handle the infrared divergences. We use the sub- 
traction method in the formulation of ref. |7I8| . 
This method treats exactly both the phase space 
and the matrix elements, thus ensuring a better 
behaved numerical result. We present a complete 
computation of jet and di-jet cross-sections and 
we compare our theoretical predictions with pre- 
vious results and with experimental data. 
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1 We use the word "infrared" both for the properly soft 
effects and for the collinear ones. 



2. Some remarks 
method 



about subtraction 



To compute a generic physical observable in 
hadronic physics it is necessary to implement the 
factorization theorem in a computer code. Thus 
we have to deal with unphysical large oscilla- 
tions (due to the cancellation of the infrared di- 
vergences following the factorization prescription) 
and the computer precision. The solution of this 
problem it is to remove divergences analitycally 
and then to use a numerical method (typically a 
MC algorithm) to integrate the finite reminders. 
In this section we briefly describe a simplified 1-D 
model [5] to outline the differences between the 
subtraction and the slicing method. 

2.1. Application of the slicing and the sub- 
traction method to a simple case 

Let us consider the integral, 



i r 1 

(F) = — dxS(l - x)F(x) 
2e Jo 



dxil-xy^Fix), 



(1) 



in the limit of e — » 0, where F is any regular func- 
tion. This is a simplified situation of what hap- 
pens to cancel out the soft emission (the second 
term in Q) and the virtual correction (the first 
term) divergent terms when we compute a generic 
infrared safe observable F; e represents the (di- 
mensional) infrared regulator. The first integral 
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is trivial, the second one is not; in general, we 
can not compute it analytically. To remove the 
e pole we need, at least, to express explicitly the 
divergent part of the second integral. If we use 
the slicing strategy we rewrite the second term of 
(JTJ as follows: 



dx- 



F(x) 



1-5 



dx- 



F(x) 



dx 



(1-x) 

F(x) 



1-5 



(1-X) 



l+2c 



< 5 < 1. (2) 



It is clear that the only divergent term is in the 
second piece on the RHS, and if we perform a 
Taylor expansion around 6 = 0, we obtain the e 
pole we need. Thus singular terms in e in Eq. 
cancel out, and we obtain 



(F) 



1 * dx Ei?l +F (i)i og 5 + o(S,e). (3) 
1 — x 



This formula is obviously finite if we set e = 
whatever 5 we choose, thus we can perform the 
numerical integration. Formally the 5 depen- 
dence of © is due only to the Taylor approxima- 
tion we have done and it is of first order because 
the logS term is canceled by an analogous factor 
coming from the upper integration limit in the 
©. Then it is clear that the smaller 6, the less 
accurate the numerical cancellations. 
On the other hand if we use the subtraction tech- 
nique we can rewrite © as 



L 



F(l)9(x-l + x c ) 



F(l) 



dx 



(1- 
6(x - 1 



(l_ x )(l+2e) 



x yi+2e) 

< x c < 1. 



(4) 



This time we have added and removed the same 
term, essentially the singular behavior of the in- 
tegrand, in such a way that the divergent piece 
is isolated in the second factor as before. But 
this time we can perform exactly the integration 
of the divergent part. After removing the e pole 
and we obtain 

f dx F(x)-F mx -l + Xc ) 



1-x 



F(l)logx c . (5) 



The logarithm term the equation (J5J is formally 
identical to the one in However, here no large 
number cancellations are present because we are 
not forced to set x c ~ 1- Thus the results from 
® are exact and free of large fluctuations, but 
a good numerical integration is still necessary to 
obtain accurate results. 

3. Structure of the computation 

As is well known, in real photon-photon scat- 
tering photons can interact directly or hadroni- 
cally. Thus the interaction is conventionally di- 
vided in three classes of Feynman diagrams: dou- 
ble resolved class, where the two photons inter- 
act by their hadronic components, the single re- 
solved class, where only one photon interact by its 
hadronic component, and the direct class, where 
the photons themselves are the initial states of the 
hard process. As we will explicitly show each con- 
tribution has no physical meaning but only the 
sum correspond to an observable. For the first 
and second class we have modified the code ^01 
used to compute photon-hadron jet observables. 
As far as direct contribution is concerned we have 
recalculated and checked with |5lllj all the rele- 
vant matrix elements, which appear in the com- 
putation at NLO (a 2 as order). Finally we have 
implemented the third class of diagrams in a com- 
puter code, analogue to the previous two codes. 
By summing the results from the three codes, we 
can obtain, using the subtraction method, any in- 
frared finite jet and di-jet observable at NLO in 
real photon collisions. The details of the compu- 
tation can be found in |12j . 

3.1. Consistency checks 

We start by presenting the comparison of our 
theoretical predictions with the LEP OPAL data 
[T3"| . In the lower panel of Figure 1 the differen- 
tial inclusive one jet cross-section is shown. As 
one can see our computation (we use in Figure 1 
MS scheme, ^ = m = = E^ ot /2, 2 GRV- 
HO photon density functions, Aqcd — 130 MeV; 
see ^21 f° r the jet definition, for the cuts and 
for the kinematic) is in decent agreement with 



2 Note that the half of total transverse energy of the event 
(E^ ot ) is equal to E T -jet at LO. 
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Figure 1. 



experimental data. In the central panel the ra- 
tios of the double resolved (dots) , single resolved 
(dashed), direct (solid) contributions and of the 
sum are plotted. Our computation agrees with 
the previous one [5] . Let us study the consistency 
of our computation, for example by changing the 
factorization scale (/i^ = 2/io, 0.5/io)- In fact if 
we change the factorization scale, upper panel of 
Figure 1, we obtain the same theoretical predic- 
tion (circles) but the double, single resolved and 
direct contributions change. Because each class of 
Feynman diagrams is separately divergent while 
the sum is finite, in order to perform the com- 
putation we had subtracted divergences analyt- 
ically from each class, adding (or removing) an 
infinite term. Such a procedure depends on the 
regularization we use, i.e. on the factorization 
scale and scheme. Hence changing regularization 
the partial results change, the total does not 3 . In 
conclusion Figure 1 represents a good consistency 
check of our computation and shows that the sep- 
aration of Feynman graphs at NLO is physically 
meaningless. 

3 Up to uncontrolled higher-order terms, which appear to 
be negligible in this case. 
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Figure 2. 



To distinguish between the two types of interac- 
tions of the photons, it is customary to introduce 
the couple of variables 

x± = Y.Jet=l,2(E J ± P j) 

where the sum in the numerator runs over the 
two hardest jets of the event. The meaning of 
the definition (JBJ it is clear: in fact if we consider 
the Born process in the direct case, we have only 
two back-to-back jets and no hadronic remnants, 
which implies x^ = 1. Thus generally the closer 
x^: to one, the more the cross section is domi- 
nated by the direct contribution. Elsewhere the 
process is dominated by resolved interactions. We 
are able to distinguish the behavior of each pho- 
ton (coming from the electron or the positron), 
and then we can study the single resolved case 
too. It is clear that we need a x^ ep value to 
separate (in an arbitrary manner) between the 
hadronic and the direct contribution: we will use 
x sep _ 0.75, the same value used by OPAL, sug- 
gested by MonteCarlo simulations pQ. Therefore 
we consider three regions on the rr^-plane, corre- 
sponding to the three types of interactions men- 
tioned before. We expect that the cross-section 
in the "hadronic" region to be the most sensitive 
to the PDF choice. Indeed this is the case. More- 
over we can predict the differential cross-sections 
with respect to x^ . To this aim it is interesting 
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to check what happens if we change again the fac- 
torization scale. In Figure 2 we plot da/dx^ for 
three different scales in the region > 0.75. For 
x± < 0.80 the curves are indistinguishable, but it 
is evident that close to one we obtain unstable 
results, strongly depending on the factorization 
scale choice. 



4. IR-Sensitive regions 

We can deduce physical observations from 
these consistency tests. In fact if we look at the 
hadronization corrections 4 taken from jlll4j to 
dajdx^j in the same region we observe that they 
are between 40% and 260%. Hence the conclusion 
is that the perturbative computation at NLO fails 
in this region. To understand this fact let us con- 
sider the relationship between x^ and the fraction 



of energy (1 



') of the emitter parton carried 



out by a soft gluon emitted in the final state of 
the direct process. It is simple to show that in a 
three-jet configuration 



„1,2 



1 



1. 



(7) 



Thus the "dangerous region" showed in Figure 2 
corresponds to a soft gluon emission kinematic, 
where the exponentiation and the resummation 
to all order of the large logarithms are neces- 
sary to obtain sensible predictions. In particular 
the dominant term of the perturbative differen- 
tial cross-section in the limit x^ — > 1 is roughly 



da 



9(x± - 1) 



1 



1 



(8) 



The expression (JSJ explicitly indicates the finite- 
ness (i.e. the infrared safeness) of the cross- 
section we expect from jet definitions, but be- 
cause infrared divergence is canceled badly (ex- 
actly for x^ = 1) large logarithms appear. Thus 
the comparison of the data in this region with 
NLO predictions cannot be satisfactory. In other 
words one has to renounce to compare the exclu- 
sive observable da/dx^ in this region with fixed 

4 Here, the hadronization correction is defined as the ratio 
between the MC cross-sections at the parton and hadron 
levels. 



order computations. To avoid this problem one 
can use more inclusive quantities |1I15| . obtained 
by integrating da/dx^ over x^ > 0.8. This im- 
proves the agreement between data and predic- 
tions. 

We also point out that there is another region in 
the phase space which can introduce large log- 
arithmic corrections to NLO computations. It 
is known that di-jet observables are sensible to 
the choice of the energy cuts on the phase space, 
in particular it is shown |10| that symmetric 
cuts over the two jets transverse energy imply 
large logarithmic corrections to certain observ- 
ables. Thus to avoid this effect OPAL collabora- 
tion has chosen asymmetric cuts on the E^ et - 
plane defined as 

E T = (#t + Et)/2 > 5 GeV, 



< 4> = 



I -^T — ^ T I 



E\ 



1 

<4- 



(9) 



Even if these cuts are free from the troubles we 
cited before, we point out that other problems 
could arise for NLO predictions. Indeed, NLO 
computation imposes cj> < 1/3, and setting the </>■ 
limit near this value could be problematic. Fur- 
thermore if the 0-limit goes to zero infrared ef- 
fects arise. In fact one should expect that the 
cross-section goes to zero if the phase space be- 
comes small, but as Figure 3 shows this is not the 
case at NLO. 
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Figure 3. Variation of the cross-section with the 
0-cut at fixed E T — 5 GeV 



The explanation of this is simple: if — ^ then 
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there is not enough phase space to emit a soft 
gluon to cancel out properly virtual divergences. 
Therefore large logarithms of <j> J2| become domi- 
nant when 0^0. Thus even if in principle OPAL 
choice £Q could be dangerous in practice it is far 
from the "dangerous regions" . 

5. Phenomenological results 
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Figure 4. Differential di-jet cross-section mea- 
sured by OPAL compared with our prediction 



dictions and data for the most inclusive observ- 
able and in the single resolved region, whereas 
the agreement does not seem satisfactory in the 
double resolved region. Even if we consider the 
theoretical errors due to the scale and the PDF 
choice or Kqcd, as shown in Figure 4, we are not 
able to adsorb discrepancy in the theoretical un- 
certainties. Furthermore the situation is worse if 
we consider da/dx^ in this region, NLO predic- 
tion cannot satisfactory describe the experimen- 
tal data. The same discrepancies affect the com- 
parison between data and the previous theoretical 
computation [Q. Thence we conclude that more 
experimental and theoretical studies are needed. 




In this section we compare our results to recent 
data. In Figure 4 our theoretical predictions are 
compared to OPAL measurements pQ (see that 
paper for the kinematics and the cuts) for di- 
jet cross-sections are plotted . The three series 
of curves correspond respectively to the whole 
x^: region, either x+ > 0.75 or x~ > 0.75, and 
x^ < 0.75. We use for the central theoretical pre- 
diction GRV-HO in DIS scheme with fio = fin = 
H F = E^ ot /2 and A QCD = 130 MeV. The NLO 
predictions, for all the observables we present, are 
divided by the factor (1 + Sh), the hadroniza- 
tion corrections taken from |H14| . These terms 
are less than 10% for E-r grater than 10 GeV, 
but rise to about 20% towards small Et- Plots 
show a very good agreement between QCD pre- 



Figure 5. Differential di-jet cross-section with re- 
spect to x± 5 < E T < 7 GeV. 



If we compare (Figure 5) da/dx^ with the data 
for < x± < 1 and 5 < E T < 7 GeV (where the 
"resolved" component is dominant) we have only 
a partial agreement. We note that the correspon- 
dence we find for x^ > 0.8 is not particularly 
meaningful. In fact, as we have discussed in Sec- 
tion 3 the result is not stable under variation of 
the factorization scale. On the other hand in the 
energy range 7 < E T < 11 GeV data agree well 
with the theoretical curves outside the "danger- 
ous region", i.e. for x^ < 0.8. Other quanti- 
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ties measured by OPAL have been computed and 
good agreement has been found with the theoret- 
ical predictions. 
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Figure 6. Differential single jet cross-section mea- 
sured by L3: Preliminary Data 



Recently, an L3 analysis has been presented |2] , 
for the inclusive single jet da/dEx on a very ex- 
tended E T range (E™ ax w 60 GeV). In Figure 6 
we present the comparison between our predic- 
tion and L3 preliminary data. As one can see it is 
possible to divide the momentum range in two re- 
gions. For momenta less than 15 GeV the theoret- 
ical curve has the right shape to fit the data but 
it has the wrong normalization (approximately it 
is a factor two higher). For large momenta data 
and theoretical prediction are clearly in disagree- 
ment. The discrepancy is at places of one order 
of magnitude, which cannot be accounted for by 
the theoretical uncertainties. In the upper panel 
of Figure 6 the scale dependence of the theoretical 
prediction (circles) is shown, and similarly to the 
case of Figure 1 the variation is very small. Dot- 
ted curve re-shows the scale dependence of the 
"hadronic" contribution. 



6. Conclusion 

We have implemented a computer code which 
allows to compute jet and di-jet observables at 
NLO using the subtraction method. Our results 
are in agreement with the previous theoretical 
predictions obtained with the slicing method 5 . 
The comparison of our results with OPAL data 1 
is satisfactory for inclusive observables, whereas 
the comparison with L3 data |2] display sizable 
disagreement in the large Et region. 
Finally, we have analyzed some theoretical dif- 
ficulties to describe particular exclusive observ- 
ables measured in pQ. The theoretical uncertain- 
ties for inclusive quantity are under control. 
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